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ABSTRACT 

This  paper  derives  the  optimal  cycle  time  and  total  cost  of  an  EPQ  model  with  a  variable  rate  of  deterioration 
and  mixed  demand  pattern.  Demand  is  inventory  dependent  during  the  inventory  buildup  time  and  constant  demand  has 
been  used  for  during  an  in  ventory  depletion  period.  The  rate  of  deterioration  is  changing  after  particular  time  period.  A 
numerical  example  has  been  provided  to  validate  the  model.  The  result  of  sensitivity  analysis  indicated  that  total  cycle  time 
and  cost  are  sensitive  to  production  rate,  increase  in  deterioration  rate  and  demand  rate. 
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1.  INTRODUCTION 

In  the  manufacturing  sector,  when  items  are  produced  internally  instead  of  being  obtained  from  an  outside 
supplier,  the  EPQ  model  is  often  employed  to  determine  the  optimal  production  lot-size  that  minimizes  overall 
production/inventory  costs.  The  complexity  of  a  resulting  EPQ  model  depends  on  the  assumptions  one  makes  about 
various  parameters  of  the  inventory  system.  Among  these  features  are  demand  types,  actions  for  excess  demand, 
production  rates,  treatment  for  imperfect  quality  items,  minimal  service  level  constraint,  various  cost  parameters, 
and  issuing  policies  when  items  are  perishable.  The  classic  EPQ  model  assumes  that  manufacturing  facility 
functions  perfectly  during  a  production  run.  However,  due  to  process  deterioration  or  other  factors,  the  generation 
of  imperfect  quality  items  is  inevitable. 

From  product  life  cycle,  demand  rate  will  be  constant  only  when  the  product  maturity  is  on  stage.  In  the 
competitive  scenario,  higher  stock  of  items  attracts  the  customers  to  purchase  more  (Teng  et.al  2005).  Deterioration 
is  the  common  phenomenon  in  the  real  production  system  of  products  like  pharmaceuticals,  food,  flowers, 
vegetables,  etc.  The  effect  of  an  imperfect  production  process,  by  assuming  different  rates  of  deterioration  during 
the  production  process,  on  the  optimal  cycle  time  had  been  studied  by  (Rosenblatt  and  Lee,  1986).  The  EOQ  model 
with  power  form  of  stock  dependent  demand  for  deteriorating  items  (Teng,  et.al  2005),  investigated  that  the  shape 
parameter  of  the  demand  and  the  selling  price  is  highly  sensitive  to  optimal  solution.  The  impact  of  the  random 
machine  failure  on  the  EPQ  model  (Teng,  et.al,  2005 )  was  that,  the  demand  function  and  purchase  cost  is  positive 
and  fluctuating  with  time.  The  production  process  can  shift  at  random  from,  in  control  to  out  of  control  state,  during 
the  production  process.  In  such  situation,  the  items  produced  in  out  of  control  state  deteriorate  at  a  higher  rate,  than 
the  normal  rate;  hence  get  consumed  by  demand  under  the  LIFO  policy  (Garry  and  Dah,  2006).  Behrouz  and 
Babak,  developed  economic  production  quantity  (EPQ)  model,  by  considering  both  the  depreciation  cost  of  stored 
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items  and  process  quality  cost.  Depreciation  cost  and  process  quality  cost  was  assumed  to  be  continuous  functions  of 
holding  time  and  of  production  run  length.  Yuan  et.al  (2011)  derived  the  optimal  manufacturing  batch  size  and  number  of 
shipments  with  scrap,  using  mathematical  modeling  and  algebraic  approach.  Ata  (2011)  discussed  EPQ  model  for  multi 
products  single  machine,  with  discrete  delivery.  Taleizadeh,  Naja  and  Naiki,  focused  on  EPQ  model  with  a  production 
capacity  limitation  and  a  random  defective  production  rate.  Gede  and  Hui  (2010),  analyzed  an  EPQ  model  for  deteriorating 
items  with  stochastic  machine  unavailability  and  price-  dependent  demand.  Jinn  (2007),  used  time  varying  demand  and 
cost,  to  analyze  EPQ  model  and  characterize  the  influences  of  both  demand  and  cost,  over  the  length  of  production  run 
time  and  the  economic  production  quantity.  Hezari  (2008)  developed  EPQ  model,  by  considering  imperfect  and  defective 
items.  The  holding  cost  for  defective  items  has  not  been  considered.  David  et.  al  (2009),  considered  partial  backordering 
with  constant  demand,  to  study  inventory  model.  Jiang  Tao  (2009)  considers  a  multi-item  inventory  system,  where  the 
vendor  provides  the  retailer  with  delay  in  payment.  Yao  (2007)  established  a  model  for  deteriorating  items  under  delay  in 
payment,  in  which  the  demand  is  a  negative  exponential  function  of  price.  Jia  (2009)  reported  a  supply  chain  system  with 
trade  credit.  Liao  [6]  established  an  EPQ  model  for  deteriorating  items  with  delay  in  payment.  Ardak  et.  al.  (2017) 
developed  EPQ  model  by  considering  mixed  demand  pattern;  inventory  dependent  demand  at  inventory  buildup  time  and 
constant  demand  during  depletion  period.  Ardak  &  Borade  (2017),  considered  time  dependent  holding  cost  to  developed 
EPQ  model. 

Most  of  the  inventory  models  considered  various  demand  patterns  like  stock  dependent  demand,  power  form 
demand,  ramp  type  of  demand,  time  dependent  demand,  selling  price  dependent  demand  and  exponential  demand. 
It  is  well  known  that,  the  demand  rate  varies  with  change  in  inventory  level.  From  the  literature  review,  it  is  observed  that 
the  different  demand  pattern  has  not  been  discussed  by  any  researcher  so  far  on  different  time  period.  The  present  study 
may  be  significant  in  filling  this  gap  since  it  aims  to  develop  EPQ  model  by  assuming  demand  as  inventory  dependent 
criteria  during  the  inventory  buildup  time  and  constant  during  an  inventory  depletion  period.  This  paper  has  five  sections. 
Research  motivation  and  literature  is  narrated  in  section  1 .  Section  2  contains  notations  and  assumption.  The  methodology 
used  to  develop  the  model  is  discussed  in  section  3.  Numerical  and  sensitivity  analysis  is  discussed  in  section  4  and  finally 
concluded  in  section  5 

2.  ASSUMPTION  AND  NOTATIONS 

Following  assumption  and  notations  are  used  to  develop  the  model. 

2.1  Assumptions 

The  following  assumptions  are  made  in  the  development  of  the  model. 

•  The  production  rate  is  known  and  is  constant. 

•  The  production  rate  is  greater  than  the  demand. 

•  The  demand  is  inventory  level  dependent  in  uptime  and  constant  in  down  time. 

•  Deterioration  of  the  items  is  varying. 

•  Inventory  holding  cost  is  known  and  termed  as  constant. 

•  Deterioration  of  the  items  starts  as  it  enters  into  inventory. 
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•  Shortages  are  not  allowed. 

•  Every  produced  item  needs  inspection. 

2.2  Notation 

Ii  -  Inventory  level  during  production  uptime,  I2  -  Inventory  level  during  production  down  time,  Ti-  Production 
up  time.  Maximum  inventory  buildup  in  this  time,  T2-Production  down  time.,  T2-Production  down  time  where  rate  of 
deterioration  changes,  P  -  Rate  of  production,  D  -  Basic  demand,  0  -  Basic  rate  of  Deterioration,  (1-  Increase  in  rate  of 
deterioration,  (-  Inventory  dependent  consumption  rate  parameter.,  H  -  Holding  cost  per  unit,  Q  _ Inspection  cost  per  unit, 
T  -  Production  cycle  time.,  TC  -  Total  cost.,  IC  -  Total  Inspection  cost,  TCT  -  Total  cost  per  unit  time,  TI  -  Total 
Inventory. 

3.  MODEL  FORMULATION 


This  EPQ  model  considered  inventory  dependent  demand  and  constant  holding  cost  and  different  rate  of 
deteriorating  items.  The  production  starts  at  time  zero  and  inventory  start  to  build  up  at  the  rate  of  P  -  D,  which  is  the 
production  rate  minus  consumption  rate  with  deterioration  up  to  time  Ti  where  the  production  process  stops  and  the 
inventory  on  hand  reaches  its  maximum  level.  AfterTi,  the  inventory  consumed  under  the  action  of  demand  and 
deterioration  up  to  time  T2.  After  time  T2,  the  inventory  level  decreases  with  the  constant  consumption  rate,  D,  but  the  rate 
of  deterioration  increases  until  it  becomes  zero  at  times.  The  objective  of  this  part  of  this  research  is  to  find  the  optimal 
inventory  cycle  time. 

Q1 

>- 
tr 
o 


Figure  1:  Inventory  System 


As  shown  in  figure.  1,  the  production  will  start  at  t  =  0,  during  the  time  period  (0,  T i)  the  inventory  will  gradually 
build  up  with  constant  deterioration.  Later,  production  stops  and  build  up  inventory  is  consumed  to  fulfill  the  demand. 
For  the  time  period  (0,  T2),  inventory  will  consumed  under  the  action  of  constant  demand  and  basic  rate  of  deterioration. 
Maximum  inventory  will  be  at  time  t  =  Tt.  For  the  time  period  (0,  T3),  inventory  gets  consumed  under  the  action  of 
constant  demand  and  increase  rate  of  deterioration.  Production  system  can  be  described  by  the  following  differential 
equations. 


dl,  (t) 

— t =  P-D-od1(t)-0I1(t)  0<t<T, 


(1) 
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During  the  time  interval  (0,  T2),  deterioration  of  items  starts  at  basic  rate.  Hence  the  system  is  affected  by  the 
combined  effect  of  constant  demand  and  deterioration.  Hence,  the  change  in  inventory  level  is  governed  by  the  following 
differential  equation. 


=  -D-0L,(t)  0  <  t  <  T2 


(2) 


In  a  time  interval  (0,  T3),  demand  rate  is  constant  and  the  rate  of  deterioration  of  items  is  at  increase  rate.  Hence, 
the  change  in  inventory  level  is  governed  by  the  following  differential  equation. 


dl3(t) 

dt 


=  -D-|3I3(t)  0 < t < T3 


(3) 


Initial  boundary  conditions  associated  with  this  equation  are,  at  t  =  0, 13  (t)  =  0,  at  t  =  T2, 12  (T2)  =  Q  and  at  t  =  T3, 13  (T3)  =0 
the  solution  for  above  equations  is  as  follows.  These  three  equations  are  used  in  the  derivation  of  our  model. 


Il(t) 


P-D 
a  +  0 


1-e 


-(a+0)t 


(4) 


(5) 


6 


By  using  initial  boundary  conditions  I3  (0)  =  Q  for  equation  6 


1 

1 _ 

^  PTf 

e  v  -1 

=  D 

(7) 


By  using  boundary  conditions  I3  (T 3)  =  I2  (0)  and  Eq.7  for  equation  4  and  5, 


P-D 


a  +  0 


-(a+0)T, 

1-e 


-D 


Q  +  - 

0 


fflPi)  „  (p-D*Ti)  T  pr 

'  2“  D  3'  2 


(8) 


Total  Inventory  is  given  by 


T|  ^2  ^3 

TI  =  J  I,  ( t ) dt  +  |  I0  (t) dt  +  j  Io(t)dt 
0  1  0  Z  0  J 


(9) 


Total  Inventory  holding  cost  is  given  by 
T  T| 

J  I[(t)dt+  J  I0(t)dt+  J  I3  (t)c 
0  1  0  z  O'5 


IH  =  H 


( 


=  H 


— xT2  +D 
2  1 


T3+ 


if 


(P-OxTi)  PT; 


2 't 


D 


-T,  - 


3"  2 


D 

+  —  xTT 
2  3 


(10) 


All  produced  items  need  an  inspection,  so  inspection  cost  is  given  by. 
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IC'Ci/V')*  =ci(^Ti2 


(11) 


From  equations  10  and  1 1 

Total  cost  =  Set  up  cost  +  Holding  cost  +  Inspection  cost. 


f 


TC  =  A  +  IH  +  IC  =  A  +  H 


^xT2  +  D 
2  1 


T3+ 


PT;r 


(P-DxTj)  PT. 


D 


-X 


3  2 


+  — xT2 
2  3 


+  C.  I  — xX.2 

n  2  i 


(12) 


Production  cycle  time  =  T  =  T|  +  T0  +  X^ 


(I3) 


(P-DxTj)  PX, 


v  , ,  ■  r3  (  P-D^i  PTf 

X  =  X,  +  - — -X, — — +x.  =X  1  +  -  1  3 

1  D  3  2  3  H 


D 


(14) 


Total  cost  per  unit  time  is  given  by 


Tcr  _  TC  _  A  +  Klf+Kffl  n[K,.Ki.K,l|]  kX  K,l 
T 


T.K,  --X? 

1  6  2  3 


(15) 


Where, 


K 


,P-D  „  TT/ _ ,  „  HD  _  HDP  HP(P-D)  HDp2  „  ,  P-D 


:(h  +  ci)^,k1=h(p-d),k2  =  ^,k3  =  ^,k 


2  ’  2  ’  ^5  4  ’  ''6 


-Ka=1  +  +- 


D 


Total  cost  per  unit  time  is  a  function  of  Ti  and  T3.  The  optimum  production  up  time  can  be  derived  from 
satisfying  the  equation. 


dTCT 

and  - - =  0  16 

dT 

Error!  Bookmark  not  defined. 

By  solving  equation  No.  15  and  16  simultaneously,  optimum  value  of  Ti  and  T3can  be  obtained. 


dTCT  d 


dT3  dT3 


^3  3  5  3 


A  +  KT2  +  KjT|T3  - T2  +  K3  +  K4T: ]  -  K,T3  -  KCT4 

T.K.  -—T^ 

1  6  2  3 


=  0 


dTCT  d 


dTj  dTj 


A  +  KTf  +  K,  T,  To  -  To  I  Ko  +  Ko  +  K ,  X  I  -  KoT3  -  KcTn4 


1  "  1  3  3  L  2 "  3  4 


1]- 


L3  3  5' 3 


T1K6-|T32 


Ignoring  higher  power  of  T3  &  Tj  and  after  some  approximation,  following  quadratic  equation  can  be  obtained. 
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AKa 

T3  =  - —~r  <17) 

VK2K6  +  K3K6-|K 

Tl2  (K1K6 ) '  T1  [2K2K6T3  +  2K3K6T3  +  KPT3] '  APT3  =  0  (18> 

Closed  form  solution  can  be  obtained  by  considering  positive  root  of  equation  18. 

4.  NUMERICAL  AND  SENSITIVITY  ANALYSIS 

Numerical  example  and  sensitivity  analysis  have  been  carried  out  to  validate  the  theoretical  aspects.  The 
numerical  data  is  adapted  from  Ardak  and  Borade.  (1)  Let,  A=  Rs.30  per  production  cycle,  P  =  2500  units  per  unit  time,  D 
=  1200  units  per  unit  time,  (=  0.5,  9  =  0.1,  a=  2  b  =  1.5.  And  P  =  0.2  the  optimum  value  of  Tt  can  be  found.  The  optimum 
value  of  U  is  1.4.  The  optimum  total  cost  per  unit  time  is  TCT  =  Rs.2411.  Sensitivity  analysis  is  carried  out  by  changing 
each  parameter  by  -40%  to  +40%,  taking  one  parameter  at  a  time  and  keeping  the  others  unchanged. 


Table  1:  Sensitivity  Analysis  of  Tj  and  TCT 


Parameter  Changes 

Parameters 

40% 

-20% 

20% 

40% 

Tj 

TCT 

T, 

TCT 

Tj 

TCT 

T! 

TCT 

P 

4.42 

1794.77 

1.99 

1986.96 

1.1 

3015.83 

0.96 

3767.80 

D 

1.06 

3745.65 

1.21 

2841.55 

1.7 

2236.12 

2.14 

2209.72 

P 

1.34 

2330.63 

1.37 

2371.09 

1.4 

2450.42 

1.46 

2489.23 

H 

1.86 

2209.04 

1.58 

2301.22 

1.3 

2525.52 

1.17 

2640.02 

Ci 

1.38 

2155.99 

1.39 

2282.64 

1.4 

2541.15 

1.42 

2673.05 

A  careful  study  of  the  computational  results  as  shown  in  Tables  5.5.6,  and  within  the  range  of  values  of  the  chosen 
parameters,  reveals  the  following  observations: 

From  Table  1,  a  higher  value  of  results  in  higher  values  of  demand  rate,  P  and  inspection  cost,  but  lower  values  of 
results  in  higher  values  of  production  rate  and  holding  cost.  This  implies  that  increase  in  the  demand  rate  will  result  in  the 
increase  in  the  length  of  inventory  buildup  time.  This  is  expected  since  if  demand  rate  is  higher,  the  stock  will  take  more 
time  to  build  up.  Increase  in  deterioration  rate  (P)  will  result  in  an  increase  in  the  optimal  length  of  time  with  positive 
inventory.  This  attracts  the  attention  of  the  inventory  managers.  With  the  deterioration  rate  higher,  however,  increases 
inventory  buildup  time.  Increase  in  production  rate  results  in  decrease  production  up  time.  During  production  up  time 
demand  used  is  inventory  dependent,  so,  as  rate  of  production  increases  inventory  increases,  but  with  an  increase  in 
inventory  demand  also  increases.  The  nature  of  curve  for  production  rate  shown  in  figure  1  indicates  this  change.  The 
effects  of  stock  dependent  demand  rate  and  deterioration  on  the  optimal  replenishment  policy  are  significant,  and  hence 
should  not  be  ignored  in  developing  inventory  models.  It  is  assumed  that  all  produced  items  needs  inspection.  Production 
up  time  is  increasing  function  of  inspection  cost.  The  inspection  cost  is  slightly  sensitive  to  production  up  time.  Figure  2 
shows  the  sensitivity  analysis  of  total  cost  per  unit  time.  Production  rate,  demand  rate  is  highly  sensitive  to  TCT.  The  rate 
of  deterioration  and  holding  cost  is  moderately  sensitive  to  the  total  cost  per  unit  time.  Increase  in  demand  reduces  the  cost 
and  increase  in  production  rate  increases  the  total  cost.  Deterioration  rate  also  increases  the  TCT. 
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Figure.2:  T,  v/s  Parameter  Change  Figure.3:  TCT  v/s  Parameter  Change 


Table  2:  Sensitivity  Analysis  of  T2  and  T3 


Parameter  Changes 

Parameters 

-40% 

-20% 

20% 

40% 

T  2 

t3 

T  2 

t3 

T  2 

t3 

T  2 

t3 

P 

0.13 

0.9 

0.17 

1.05 

0.2 

1.3 

0.23 

1.41 

D 

0.23 

2 

0.22 

1.49 

0.2 

0.98 

0.14 

0.83 

P 

0.16 

1.21 

0.18 

1.2 

0.2 

1.17 

0.23 

1.16 

H 

0.23 

1.54 

0.21 

1.33 

0.2 

1.08 

0.17 

1.08 

Ci 

0.18 

1.18 

0.19 

1.18 

0.2 

1.19 

0.21 

1.19 

Figure  4:  T2  v/s  Parameter  Change  Figure  5:  T2  v/s  Parameter  Change 


Values  from  table  2  and  Figure  4  shows  the  sensitivity  analysis  of  production  down  time  T2.  Time  T2  is  highly 
sensitive  to  production  rate,  demand,  holding  costs  and  rate  of  deterioration  ‘P’.  An  increase  in  the  demand  rate  will  result 
in  a  decrease  in  the  length  of  downtime.  This  is  expected  since  if  demand  rate  is  higher,  the  stock  will  take  less  time  to 
finish.  Increase  in  value  of  the  rate  of  deterioration  ‘P’  increases  the  inventory  depletion  time.  With  the  deterioration  rate 
higher,  the  stock  will  take  more  time  to  finish.  Increase  in  deterioration  rate  follows  linear  relation  with  the  increase  in 
production  time  T2.  Though  the  rate  of  deterioration  is  basic  during  T2,  but  effect  of  the  increase  rate  of  deterioration  is 
dominating.  Increase  in  rate  of  production  increases  the  inventory  and  hence  required  more  time  to  finish.  The  inspection 
cost  is  also  increasing  function  of  down  time. 


Table  2,  shows  the  sensitivity  analysis  of  T3.  This  time  period  is  highly  sensitive  to  demand  rate,  production  rate 
and  holding  cost.  During  production  down  time  demand  used  is  a  constant  demand.  Higher  values  of  constant  demand 
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gives  lower  values  of  production  downtime.  Increase  in  demand,  decreases  the  down  time  and  hence,  optimum  cycle  time 
also.  During  T3  rate  of  deterioration  is  assumed  to  be  increased.  Figure  5  indicates  that,  this  increase  rate  of  deterioration 
decreases  the  down  time.  The  increased  rate  of  deterioration  reduces  the  inventory  of  good  items  and  hence  the  time 
required  to  deplete  this  good  item  inventory  decrease.  This  may  led  to  loss  to  organization.  As  usual  increased  rate  of 
production  increases  the  time  required  to  deplete  build  up  inventory.  The  inspection  cost  is  slightly  sensitive  to  T3. 

CONCLUSIONS 

In  this  study,  the  theoretical  EPQ  model  has  been  developed  over  an  infinite  time  horizon.  Some  interesting 
observations  are  presented.  During  the  inventory  buildup  time,  demand  used  was  inventory  dependent  and  during  an 
inventory  depletion  time  it  was  constant  Sensitivity  analysis  shows  that  TCT  as  well  as  production  up  time  are  highly 
sensitive  to  production  and  demand  rate.  Increase  in  rate  of  deteriorations  increases  the  TCT  and  optimum  cycle  time.  This 
model  can  be  useful  for  the  inventory  managers  in  decision  making,  especially  for  the  perishable  items.  The  model  can  be 
further  developed  by  considering  different  deterioration  rate,  production  rate. 
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